Signatures of granular microstructure in dense shear flows 
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Granular materials react to shear stresses dif- 
ferently than do ordinary fluids. Rather than 
deforming uniformly, materials such as dry sand 
or cohesionless powders develop shear bands: 
narrow zones containing large relative parti- 
cle motion leaving adjacent regions essentially 
rigid Since shear bands mark areas of 

flow, material failure and energy dissipation, they 
play a crucial role for many industrial, civil engi- 
neering and geophysical processes [^. They also 
appear in related contexts, such as in lubricat- 
ing fluids confined to ultra-thin molecular lay- 
ers Detailed information on motion within a 
shear band in a three-dimensional geometry, in- 
cluding the degree of particle rotation and inter- 
particle slip, is lacking. Similarly, only little is 
known about how properties of the individual 
grains - their microstructure - afftect movement in 
densely packed material j^. Combining magnetic 
resonance imaging, x-ray tomography, and high- 
speed video particle tracking, we obtain the local 
steady-state particle velocity, rotation and pack- 
ing density for shear flow in a three-dimensional 
Couette geometry. We find that key characteris- 
tics of the granular microstructure determine the 
shape of the velocity profile. 

In order to probe the role of microstructure inside the 
narrow granular shear zone, independent determinations 
of the velocity and density profiles with spatial resolu- 
tion well below the size of individual particles are re- 
quired. Non-invasive measurements of this type so far 
have been limited to two-dimensional (2D) geometries 
where optical tracking of all particle positions is straight- 
forward Pt^||^-[Tll| ■ In a 3D Couette cell, as sketched in 
Fig. la, a steady-state shear flow can be set up by con- 
fining granular material between two concentric, vertical 
cylinders and turning the inner cylinder at constant ve- 
locity Uwaii while keeping the outer wall at rest. Unlike 
2D Couette cells [P|JlO||, where particles are confined to a 
single layer with constant volume, there is a free upper 
surface allowing the packing density to adjust via feed- 
back between shear-induced dilation and gravity. 

The difficulty of imaging the interior has restricted 
studies of 3D granular Couette systems to either probing 
only the surface to tracking colored tracers in very 
narrow (few particles wide) cells p^ ], or to measuring 



global quantities, such as the total applied torque and its 
fluctuations . Here we use magnetic resonance imag- 
ing (MRI) to obtain flow velocities from the interior of 
a 3D system We have used oil-rich seeds as a 

source of free protons that can be traced using MRI . 
Two kinds of seeds were used to explore the role of mi- 
crostructure: mustard seeds (spherical with mean diam- 
eter d = 1.8mm) and poppy seeds (kidney-shaped with 
mean diameter d = 0.8mm). The wall friction was con- 
trolled by gluing a layer of seeds on both cylinders. Us- 
ing a spin-tagging technique, horizontal slices were im- 
aged, as sketched in Fig. la. In the resulting MRI image 
(Fig. lb), the shear band shows up as the narrow re- 
gion of deformed stripes near the inner, moving, cylinder. 
Imaging slices at different heights, h, we measured these 
deformations from which the azimuthal velocity profiles 
were calculated throughout the cell. Similarly, three- 
dimensional x-ray tomography allowed us to calculate 
packing fraction profiles at various heights (Fig. Id). At 
the transparent cell bottom additional velocity and pack- 
ing fraction information was gathered by direct, high- 
speed-video particle tracking (Fig. Ic). 

Prior to each set of measurements, the cell was run 
until steady-state was reached (as determined by the 
fact that the packing density profiles became station- 
ary). Over a total acquisition time of typically ITmin 
per MRI slice, high-resolution images revealed the long- 
time average local displacements during a time interval, 
St = 100 ms (Fig. lb). The azimuthal velocity at a given 
distance r from the inner wall was obtained by exploit- 
ing the cylindrical symmetry of the Couette geometry: 
we calculated the MRI intensities along a circle of radius 
r, /turning (6) and /stopped (6), for the cell turning and at 
rest, respectively. The position of the central peak in the 
cross-correlation between /turning(0) and /stoppcd(0) cor- 
responds to the average azimuthal distance traveled by 
the material during St. This technique yielded the angle- 
and time-averaged radial profile of the steady-state az- 
imuthal mass flow velocity, v(r), resolving v{r) to within 
O.lmm/s and r to within 0.1mm. Note that this high- 
resolution mass flow velocity is the average velocity of 
all material at radius r and thus not only contains in- 
formation about particle translation and spin, but also 
depends on packing fraction. This differs from the aver- 
age velocity of particle centers that typically is obtained 
by video particle tracking techniques. 
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FIG. 1. Non-invasive MRI, x-ray tomography and 
high-speed video probes of granular Couette flow. (a) 
Sketch of the Couette-type shear cell consisting of two 
concentric cylinders with diameters 51mm and 82mm and 
filled to a level of 60mm. The inner cylinder was ro- 
tated at angular velocities from 0.6-45rpm, corresponding to 
0.6mm/s < Wwaii < 120mm/s. The flow velocity was mea- 
sured using a MRI spin-tagging technique jlj,^. Before 
imaging, proton spins were encoded (spin-tagged) so as to dis- 
play parallel stripes when imaged. Images were taken of 5mm 
thick horizontal slices at various heights, h. (b) When the 
inner cylinder was rotated, distortion of the stripe pattern re- 
vealed the displacement of the material that occurred during 
the 100ms interval between spin-tagging and imaging. 2048 
spin-tag-image steps were used to assemble each complete im- 
age, (c) High-speed video frame, taken at 1/lOOOs, of mustard 
seeds observed through the cell's transparent bottom. To indi- 
cate the movement history of individual particles, their center 
positions over the preceding 200 frames are traced by black 
lines. The particle coloring reflects the magnitude of each par- 
ticle's average velocity during this 0.2s interval: fast particles 
(yellow) near the inner wall appear to move smoothly, while 
slower particles (orange and red) display more irregular and 
intermittent motion. In the long-time average, measured by 
the MRI technique, the flow appears smooth everywhere, (d) 
Two-dimensional, horizontal slice through the cell at half the 
filling height, taken from x-ray tomography data set. Par- 
ticles, in this case mustard seeds, appear bright (as do the 
walls of the cell) . These measurements were performed at the 
Advanced Photon Source at Argonne National Lab in a Cou- 
ette cell of slightly smaller size but under otherwise identical 
shearing conditions. 



We found v{r) highly reproducible from run to run 
even though on shorter time scales there are rapid veloc- 
ity fluctuations from point to point within the material 
that, at the boundaries, can be observed with high-speed 
video (Fig. Ic) jl^. Within the resolution of our mea- 
surements, v{r) for both types of seeds did not vary with 
height, including the regions near the top and bottom 
surfaces (Note that this indicates that non-azimuthal, 
secondary flow was small and did not affect v{r), in ac- 
cordance with tracer-bead studies jl^). Aside from an 
overall scale factor, we also did not detect any shear-rate 
dependence to the velocity proflle over the entire range 
5mm/s < Wwaii < 120mm/s explored with MRI. (With 
video imaging we verifled rate independence down to ve- 
locities as small as 0.6 mm/s) . 

For (nearly) monodisperse smooth, spherical particles, 
the decay of v(r) away from the shearing wall is dom- 
inated by abrupt drops at integer mifltiples of r/d as 
shown in Fig. 2a for mustard seeds. In the normalized ve- 
locity gradient 7(r) = (Fig. 2b) these drops appear 
as deep narrow valleys. They are correlated with pro- 
nounced oscillations in the packing density, p{r), which 
signal the presence of well-defined, single-grain- wide lay- 
ers near the moving wall (Fig. 2c). (We note that the av- 
erage density approaches the random close packing value 
at larger r in a manner that is consistent with an expo- 
nential form ]lO|.) The shear-induced layering is remi- 
niscent of that seen [|l^ in the coUisional regime of dilute 
flows, but occurs here in the high-density limit of rate- 
independent, frictional flow. Since the velocity drops are 
highly localized, we associate them with slipping at the 
interface between adjacent layers. Either directly from 
v{r) or by comparing the integrated areas of each valley 
in 7(r) we find that across each slip zone the velocity de- 
creases by approximately the same factor, b = 0.36±0.13. 

In addition to slip between layers, MRI resolves a non- 
zero velocity gradient j{r) within each layer. This can be 
caused by either particle rotation or disorder in the layer- 
ing along the radial direction. (Along the azimuthal and 
axial directions, particles within layers certainly show 
packing disorder as can be seen from the maxima in p(r) 
in Figs. 2c and 3c whose values lie significantly below 
those for a crystalline configuration). For perfectly ar- 
ranged layers, 7(r) at the layer centers would be deter- 
mined solely by particle rotation (spin) within each layer. 
The presence of particles which do not lie perfectly within 
the layers (revealed by the non-zero values of p{r) be- 
tween layers in Fig. 2c and seen directly from the tracks 
in Fig. Ic) reduces the gradient in the slip region and 
increases it at the layer centers. In the limit of complete 
disorder, we would expect the staircase shape of v{r) to 
vanish completely. Using information about the density, 
Pcir), and velocity, Vc{r), of the particle centers, as de- 
termined by the high-speed video and x-ray experiments, 
we find that approximately 90% of 7(r) at the layer cen- 
ters is due to the radial disorder. In order for the MRI 
and high-speed video experiments to give the same 7(r), 
particle spin must be considered, revealing the spin pro- 
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file seen in the inset of Fig. 2a. We note that the spin is 
small, increasing slowly with distance from the shearing 
wall to a value at r/d = 4.5 of less than one full particle 
rotation per 13d translation. 




0.4 



1 2 3 4 5 6 

r/d 

FIG. 2. Radial velocity, spin and packing density profiles 
for spherical mustard seeds, (a) Steady-state, angle- averaged 
azimuthal velocity v{r) across the shear band at h = 30mm, 
halfway below the filling level. The layer of seeds glued to 
the inner cylinder wall extends to the dotted vertical line 
at r/d = 1. The lower inset shows a photograph of the 
seeds. The upper inset shows the normalized particle spin 
rate, uid/v, as a function of distance from the moving wall 
(where u is the angular velocity of a particle about its center) . 
(b) Normalized velocity gradient 7 = {d/v)dv/dr for the data 
in the main panel of (a), (c) Angle- averaged, steady-state 
radial packing density profile p(r) computed from x-ray to- 
mography data as in Fig. Id, measuring the volume fraction 
occupied by seed material. For < r/d < 1 the seeds glued 
to the wall contribute to p. 

The message of Fig. 2 is that the overall shape of 
v{r) across the shear band can be understood as aris- 
ing from two main contributions to the velocity gradi- 
ent 7(r): a slip contribution in the presence of layering, 
and a second contribution associated with radial disorder. 
These two pieces are distinguished by the characteristic 
r-dependence they produce in v{r) on length scales larger 



than a single grain. A constant slipping fraction, 
7siip('')'^'' = -b, 

interface 

across each interface, by itself, leads to a velocity profile 
with an exponential decay, v{r) = vq exp{—hr / d) . The 
constant inter-layer slip appears on a background due to 
radial disorder, 7disorder(?') — — 2c(r — ro)/d, that starts 
at rg within the glued-on layer and increases linearly in 
strength with slope 2c as indicated by the dotted line 
(Fig. 2b). Such linear r-dependence in 7(r) results in a 
Gaussian profile, v{r) = voexp[—c{r / d — ro/d)^]. The 
sum of the two contributions to 7(r), after integration, 
leads to a product of an exponential and Gaussian term 
for the overall, averaged velocity profile: 



v{r) — voexp 



. Fits of this function to mustard seed data yield b = 
0.36 ± 0.13, c = 0.06 ± 0.03 and ro/d = 0.6 ± 0.8 inde- 
pendent of height and shear rate. 

The decomposition of 7(r) allows the quantitative 
tracking of the slip inside the shear band for different 
particle types. We find a greatly reduced slip rate b 
when smooth spherical particles are replaced by rough- 
ened ones ll^. More dramatic differences result from 
changes in the particle shape. For kidney-shaped poppy 
seeds we observe a comparatively smooth overall ve- 
locity profile (Fig. 3a). Much smaller modulations in 
7(r)(Fig. 3b) together with little layering in p{r) (Fig. 3c) 
indicate that interlayer slip is greatly reduced. The lin- 
ear trend 7(r) — — 2c(r/d — ro/d) is still present, how- 
ever, and leads to the Gaussian shape for v{r) seen in 
Fig. 4 over several orders of magnitude. These data also 
demonstrate explicitly that the shear rate fixed by fwaii, 
while setting the overall scale, leaves the shape of the 
profile unchanged. This allows for a collapse of all poppy 
seed data for different shear rates and heights in the cell 
onto a single Gaussian described by c = 0.11 ± 0.02 and 
ro/d = — 0.1±0.5. Layering can also be suppressed and 
radial disorder be promoted in sphere packings if wide 
particle size distributions are used. We find that these, 
too, produce essentially Gaussian velocity profiles [ p^ . 

From these results, the Gaussian component of the ve- 
locity profile emerges as a robust, generic feature; slip, 
if it occurs, is seen to augment v{r) by an additional, 
exponential factor. In contrast to the translationally in- 
variant slip, the explicit dependence of the Gaussian on 
the distance from the shearing wall indicates correlations. 
There are several possible mechanisms for radial corre- 
lations JPiP:^ 18-2^, including the formation of stress 
chains ||^ |^ ^ or particle clusters ||l8|,0 , both of which 
require high packing densities and the absence of slip or 
some degree of interlocking. However, a detailed mech- 
anism that would explain the observed Gaussian depen- 
dence presently does not exist. A successful theory must 
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FIG. 3. Radial velocity, spin and packing density profiles 
for aspherical poppy seeds, (a) Steady-state, angle- averaged 
azimuthal velocity v{r) across the shear band a,t h = 30mm. 
The inset shows a photograph of the seeds, (b) Normalized 
velocity gradient for the data in (a), (c) Corresponding an- 
gle-averaged, radial packing density profile p{r) from x-ray 
data. For < r/d < 1 the seeds glued to the wall contribute 
to p. 
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FIG. 4. Rate-independence of the velocity profiles. Log 
plot of poppy seed velocity profiles v{r) as a function of (r/d)^ 
for four different shearing wall speeds (Uwaii = 120, 52, 33, 
and 14mm/s, top to bottom). Straight lines on this plot cor- 
respond to a Gaussian with origin at r = 0. The slope, c, of 
these lines corresponds to half the average slope of the veloc- 
ity gradient in Fig. 3b and characterizes the width a = d/^/c 
of the shear band. 



slow shearing rate regime there is a direct connection be- 
tween the microstructure and the shape of the velocity 
profile. For equal-sized spherical particles, where con- 
siderable layering is found, a strong exponential contri- 
bution is observed. For aspherical seeds, which show no 
pronounced layering, the profile is almost completely de- 
scribed by a Gaussian centered on the shearing wall. Re- 
markably, despite the complicated nature of the grain- 
grain interactions, the steady-state behavior exhibits ro- 
bust behavior: the shape of the velocity profiles charac- 
terized by the two length scales, A = d/b and a ~ d/y/c, 
is found to be height and shear-rate independent. 



not only account for the linear increase in the magnitude 
of the velocity gradient with r, but also show how the 
Gaussian shape is independent of the presence of layer- 
ing. We speculate that both of these might be provided 
by a coupling of 7(r) to the average packing fraction 
and its gradients (for 2-D Couette systems a coupling 
to the average density alone has been suggested |p|,p^). 
A Gaussian, then, arises naturally for a coupling of type 
7(r) (X p{rQ) — p(r), considering first order terms in r. Be- 
cause all stress-loading is driven by the seed layer glued 
to the inner wall, its location, tq, serves as a spatial ref- 
erence point. This argument carries over to the absence 
of a Gaussian contribution to v{r) in less dense types 
of flow, where interlocking may be less effective, such 
as rapid flows down inclines with free upper surface for 
which power law forms have been reported or 
vertical gravity-driven flows through 2D "pipes" , which 
appear to follow essentially exponential profiles . 
Our results show that in this high packing density and 
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